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1. Introduction 

As it is well-known, in the improvement approach a la Symanzik [|]] the lattice QCD action has 
to be provided with an extra irrelevant contribution, the so-called Sheikholeslami-Wohlert term p]. 
In Perturbation Theory (PT), it features a scalar coefficient csw which can be Taylor-expanded in 
powers of the bare coupling go as 

csw = 4w + cswd + c §So + 0{st) • (1-D 
The zero- and one-loop coefficients have already been determined for different lattice actions [O] [ffl] 

(2) 

while c sw is still unknown: the final aim of this project is precisely to estimate it by combining the 
Schrodinger Functional formalism (SF) and the PCAC relations in the same spirit as and [||] 
where c^L and c jl were successfully recovered. 

The main difference with these two latter seminal papers lies in the fact that observables are 
evaluated perturbatively without following a diagrammatic approach but rather by means of Nume- 
rical Stochastic Perturbation Theory (NSPT), a computer algorithm characterized by a Langevin- 
like evolution of the system. 

2. Theoretical aspects - part I (basics) 

The lattice formulation of QCD we adopt is that of Wilson: a concrete expression of the well- 
known contributions to the action - namely the gauge (So), fermionic (Sp ) and Sheikholeslami- 
Wohlert (Ssw) term - can be found in [||] whose notations and conventions inspire nearly all the 
formulae appearing in this and the next section 1 . 

(2) 

A suitable observable to study in order to evaluate c sw is provided by the quark mass m q which 
can be conveniently computed by means of the lattice PCAC relation reading, 2 

l -{d« + d L )(A »O) = 2m q {P b (n)0) , (2.1) 

where O is any product of fields located at nonzero distance from n, dfi (9q) is the lattice right 
(left) derivative in the time direction and 

N f 1 N f . 

Ain) = £ V^te.^W > p \n) = tv f (n)Y5^ b fg ¥ 8 (n) , (2.2) 

f,g f,g 

where r h is a matrix acting on flavour degrees of freedom 3 . 

(2) 

In order to fix c s ^, one requires m q to be independent of contributions of order a: however, to 
achieve full improvement Eq.(2.1) has to be modified to, 

\{d^ + dlf){Al{n)Q)+c A d^d§{P h {n)0) = 2m q (P b (n)0) , (2.3) 



More generally, we stick to the setup outlined in sections 2, 4 and 6 of [Q]. 

2 From now on, the time direction will be assigned the subscript 0. 

3 Spin and colour subscripts will be usually left implicit in order to ease the notation. 
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where c a is a second improvement coefficient which, just like csw> can aiSO be decomposed as 
ca = + c^gQ + c^gQ + ^(go)- Once again, the first unknown contribution is at two-loop 
level: see [j|] and [^] for the determination of c A °^ and c A . 

The second main theoretical ingredient of the present strategy is given by the Schrodinger 
Functional: assuming the time coordinate ranges from to T and labelling the space coordinates 
as n, it consists of replacing the usual periodic boundaries by Dirichlet conditions along the time 
direction, namely, 

f*(»)Lo=o- W k (n) , Uk (n)\ no=T ^ Wl(n) (£ = 1,2,3) , (2.4) 
for the gauge degrees of freedom 4 and (P± = (I ± 7o)/2 with I being the identity matrix) 

V'WUo - P 7 (") = P + V f (n) L 0=0 , V'W |„ 0=r - P ;/ («) = P+ ¥ f (n) \ nQ=T , (2-5) 
V%)L =o "> P^") = /, + V 7/ («)L =o . V 7/ (")L =r ^ ^) =p +V ?/ (")L =r > ( 2 - 6 > 

for fermions: boundary fields W, W', p,~p, p' and p' will be defined later on. 

Due to the Schrodinger Functional formalism, the three contributions to the lattice QCD action 
get modified as follows: 

• the gauge part So becomes 

S G = /3£ nw{n)(\-^[Upy{n) + Ul y (nj\) , (2.7) 

n,fi,v V / 

H>v 

where the weight <0^ v (n) for the lattice plaquette U^ v (n) is 1 everywhere except for the 
spatial plaquette at n$ = and n$ = T whose W jlIV («) reads ^; 

• the fermionic part Sf remains in principle unchanged; anyway, in order to have one more 
parameter to play with, an additional phase e lQ v/ L n is introduced in the definition of the 
lattice covariant derivatives within the Wilson-Dirac operator: in practice, gauge fields U^(n) 
appearing in Sp are replaced by, 

U^(n) ^ e ie ^U^n) , (2.8) 

with O = and -% < d k < K for k = 1,2,3; 

• the clover term is set to for all those lattice points with hq = or no = T. 

4 Gauge fields along the time direction, defined for < «o < T, have no constraints on them. It turns out that W and 
W' can sloppily be written as W = £Pe^ c and W = £?e^ c - see section 6 of Q for notations and a more careful and 
detailed treatment of this topic - where C and C' play a similar role as the background field in classical physics: in what 
follows we will refer to the case C = C = as the trivial background. 
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3. Theoretical aspects - part II (details) 



(2) 

Before outlining the procedure that should lead to an estimate of c^, let us give a precise 
shape to the observable O appearing in Eq.(2.3): a convenient choice reads, 



Nj 

I 

f,g m,m' 



(3-1) 



where 



ft 



m 



8p f {m) 



ft 



m 



8pf(m) 



(3.2) 



After first plugging Eq.(3.1) into Eq.(2.3), then letting the derivatives with respect to p and 
p act on the Boltzmann factor and finally setting all the fermionic boundary fields to zero, some 
algebra allows one to write 



with ' 



m a 



[M fp] 



/a 12 ( H [(m+6) coc, nee] (W))eP *fg ( P -) W(J ^+6) ac, nfie] T « ^ G 
m,m' 

f D = — V (M lf T h (P \ J gk T h \ 

JP 12 ^ ^ [{m+Q)ac,nee] fs \' ~)wcj J [(m'+6) ac, nee] hl 'o ' 

tri ill' 



(3.3) 



(3.4) 
(3.5) 



with 



[(m+0) w c , n e e] 



[(m'+0)oc,nPe] 



U (m) 
U (m f ] 



cb 



- cd 



AT 



M 



-l* 



if 

[(m+6) cob, nee] 
gh 

[(m'+6)od,nPe] 



(3.6) 
(3.7) 



where M is the overall fermionic opearator in the lattice action. 

/a, fp and m q depend on the lattice spacing a, the lattice extents L^, the bare coupling go, the 
gauge fields W and W', the angles 6k (from now on, we will set the latter equal to a common value 
6) and the improvement coefficients: recalling that the approach is perturbative, we can write 6 , 



m q (L,d,x ,go,a) =m q Q \L,d,x ,a)+m q z '(L,d,xo,a)gQ + m [ g'(L,d,xo,a)gQ + ff(gQ) , (3.8) 



5 The subscript "G" stands for the mean over gauge degrees of freedom. Here and in Eqs.(3.6)-(3.7) repeated indices 
are summed over. Moreover, from now on we tacitly assume that all quantities are rescaled with a to be dimensionless. 

6 We make the dependence on W, W' , c$w and ca implicit not to overwhelm the notation; at the same time, we drop 
the subscript on the lattice extents for a reason that will become clear soon. 
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and in turn, thanks to dimensional analysis 

^M,*,*)^^^ • (3-9) 

L, Xq L, 

This formula can actually be simplified by setting the L^'s to the same value L, putting Lo = 2L 

(k) 

and choosing no = L: thus, the corrections in a to m q ' will be grouped together into a single one 
proportional to a/L. Since the aim of improvement is to get rid of lattice artifacts of order a, it 
is reasonable to estimate by requiring the only coefficient d(c^ k \c A ~ k ') left in the formula 

(2) 

above - after its reduction - to vanish. This can be achieved by the following steps: 1) fix and 
cy arbitrarily after setting c[2L c£l, c^ and to their known values; 2) perform simulations 
for different lattice extents keeping 6, W and W' constant; 3) fit the coefficient d(c^,cy); 4) 

(2) (2) 

repeat the previous steps for different choices of and c A , 5) collect the various estimates of 
d(c^,c A ^) and interpolate the values of and c A ^ for which d(c^,c^) vanishes. 

Before ending this section, some remarks are in order. 

The first term on the r.h.s.of Eq.(3.8) should normally correspond to the bare mass Mo appea- 
ring in Sp ; however, in the present setup, this is the case only if 6 = 0: we chose to set Mo = but 
to work with non- vanishing 6 to avoid any infrared divergence. 

Second, in Eq.(3.9) it is understood that mass counterterms - depending on c$w EfTl] - are sub- 

(k) 

tracted. Otherwise m q would not be in the large L limit: this subtraction prevents extra improve- 
ment coefficients to appear (see section 3 in [Q]) but, in practice, this should really matter only 
when working with renormalized quantities (while we deal with their bare counterparts). 

Finally, it is possible to disentangle the effects of a™, and by means of W and W': in 

(2) 

particular it turns out that, if the trivial background (see footnote 4) is set, only c A has an effect 
at two-loop level. We start with this choice of the boundary gauge fields to fix this coefficient, 
afterwards W and W' will be changed to determine cQ, thanks also to the by-then-known estimate 
of tf. 

4. Numerical aspects 

Two more issues have still to be addressed about the present strategy, namely how configura- 
tions are generated and how the Wilson-Dirac operator is inverted to compute /a and fp eventually: 
to answer both, we must introduce some basics of NSPT 7 . 

Its core is given by the Langevin evolution equation that, for lattice gauge variables 8 , reads 

^U^n,t) = -iY,T A [y n ^^ A S[U)+^{n,t)]U^n,t) , (4.1) 
at A 

where t is an extra degree of freedom (which can be thought as a stochastic time), S is the part of the 
lattice action depending on the U's, r\ is a Gaussian noise while V stands for the group derivative 



7 See [g| and references therein for more details on this section in general. 

8 As usual, fermion fields are integrated out so that only gauge degrees of freedom have to be eventually treated. 
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defined as (index "A" is summed over), 

^[e iaATA U^(n),U'} = ^[U^n),U'}+a A V n ,^ A ^[U^n),U'] + ... , (4.2) 

where T A are the generators of the algebra and & is a generic scalar function of both the variable 
UJn) and some more labelled U' for short. 
Given this setup, it can be shown that 

Z- l J [DU]0[U]e- s W = \imj jUt {0[U^))\ , (4.3) 

where Z is the partition function and 0[U] a generic observable depending on the gauge fields. 

Perturbation theory enters into play by formally expanding each gauge degree of freedom in 
powers of fi^ 1 - defined as /3o = 2N c /gl being N c the number of colours - up to a given order s as 

[/ /i («,f)=I+£j3 " l f/f M, (4.4) 

and then plugging this Taylor series 9 into Eq.(4.1): this results in a consistent hierarchical system 
of differential equations which can be numerically integrated by discretizing the stochastic time as 
t = m with n integer. In practice, the system starts from an arbitrary configuration and evolves 
by means of the solution of the discretized counterpart of Eq.(4. 1): the desired observable is then 
obtained by averaging its measurements on its plateau - recall the limit in t in Eq.(4.3) 10 . 

As for the inverse of the fermionic operator, the entries needed to get /a and fp can be com- 
puted by means of the following perturbative formulae 

P' 1 , 



where only the zeroth order of M has to be truly inverted: its expression for trivial W and W' can 
be found in section 3.1 of [Q]. 

5. Preliminary results 

To test the correctness of the overall setup, we computed the one-loop contribution to m q 
without any counterterm subtraction for different choices of 6 and c™- 1 1 and compared the results 

9 Strictly speaking, Eq.(4.3) is valid only if the boundary gauge fields are set to the identity as in this first part of the 
study; once that a non-trivial background field is introduced, the expansion would read Uu(n,t) = exp[(C' k —C^/T] ■ 

■ [i + Y.kP^ 1 U^ {n,t)] - consult section 6.2 in [||] for the meaning of the first term in this product. 

10 This relation is true only for continuous t so that simulations with different T values have to be performed in order 
to extrapolate to T — > afterwards. 

11 This is indeed the only c$w contribution that enters into play at this order with trivial W and W'. 
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with the analytical values in Table 1 . 
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c (0) - 


c (0) - 1 


c (0) - 1 5 


1.40 


2.67621(4) 


1.67151(2) 


0.94999(1) 


1.00 


2.63837(3) 


1.64808(1) 


0.93229(1) 


0.45 


2.60727(3) 


1.62694(1) 


0.91948(1) 


0.00 


2.60571 


1.62045 


0.91067 



Table 1: Numerical results for m^ on a 10 3 * 21 lattice with c^ = cp — 0: the last hne contains the 
infinite-volume results obtained from [^[ . 

It is reassuring that, when varying c^-, outputs change accordingly: the still-existing gap is 
explained by recalling that finite-size effects are still present and that the analytical results corre- 
spond to ntg = while in our simulations rriq ^ due to the non-vanishing values of 8 (piq 1 
approaches with decreasing 6 12 the analytical infinite-volume values computed with 8 = 0.0). 

6. Conclusions and acknowledgements 

According to the first, preliminary results, the outlined approach seems to be feasible: how- 

(2) (2) 

ever, since different extrapolations (in z and L) and interpolations (in c A and cLp when dealing 
with non-trivial W and W') are needed, extra care will have to be paid not to spoil accuracy. 

We warmly thank LRZ centre (Munich) and ECT* (Trento) for providing us with computer 
time on their clusters. 

This work was supported by the DFG SFB/TR 55. 
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